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Abstract: The distribution of binding agents (i.e., soil organic carbon (SOC) and glomalin-related soil 
protein (GRSP)) in soil aggregates was influenced by many factors, such as plant characteristics and soil 
properties. However, how these factors affect binding agents and soil aggregate stability along a climatic 
gradient remained unclear. We selected the Robinia pseudoacacia L. forests from semi-arid to semi-humid of 
the Loess Plateau, China to analyze the plant biomass, soil physical-chemical properties, SOC and GRSP 
distribution in different sized soil aggregates. We found that from semi-arid to semi-humid forests: (1) the 
proportion of macro-aggregates (>0.250 mm) significantly increased (P<0.05), whereas those of 
micro-ageregates (0.250—0.053 mm) and fine materials (<0.053 mm) decreased and soil aggregate stability 
was increased; (2) the contents of SOC and GRSP in macro-agegregates and micro-ageregates significantly 
increased, and those in fine materials decreased; (3) the contribution of SOC to soil aggregate stability was 
greater than those of total GRSP and easily extractable GRSP; (4) soil properties had greater influence on 
binding agents than plant biomass; and (5) soil aggregate stability was enhanced by increasing the contents 
of SOC and GRSP in macro-ageregates and soil property was the important part during this process. 
Climate change from semi-arid to semi-humid forests is important factor for soil structure formation 
because of its positive effect on soil aggregates. 


Keywords: binding agents; glomalin; organic carbon; soil property; soil aggregate stability 


1 Introduction 


Most ecological processes (e.g., soil respiration, gas emission and nutrient accumulation) occur in 
soil structure composed of aggregates, and the soil stability depends on the type and content of 
binding agents. Climatic conditions usually change the binding agents by affecting various 
ecological factors, such as plant characteristics and soil properties, and alter soil aggregate 
stability (Zhang et al., 2012). 

Studies have found that soil organic carbon (SOC) and glomalin-related soil protein (GRSP) 
are binding agents in soil aggregates. However, which one is more important in the formation of 
aggregates is currently controversial. Most studies believe that SOC has higher aggregation 
efficiency than GRSP, because of its higher content (Six et al., 2002). For example, Xiao et al. 
(2020) found that soil aggregate stability was strongly related with SOC along a gradient of 
secondary succession, and SOC content was higher than GRSP. Studies found that SOC was very 
inefficient in the formation of aggregates, because SOC contains several inert components (Zhang 
et al., 2012; Singh et al., 2016). On the contrary, other researchers have revealed that GRSP, 
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especially the easily extractable GRSP (EE-GRSP), is highly active, and the ability of 
glomalin-binding aggregates may be greater than that of SOC (Spohn and Giani, 2010). EE-GRSP 
in black soil of Northeast China is responsible for stabilizing soil aggregates instead of SOC 
(Zhang et al., 2012). Besides, a strong relationship between EE-GRSP and aggregate stability was 
found among different conservation tillage systems, which also supports the important role of 
EE-GRSP in the formation of aggregates. Recent studies have found that GRSP in different sized 
aggregates has varied contribution to aggregate stability (Zhang et al., 2014). For example, 
non-tillage management promotes the production of glomalin in macro- and micro-aggregates, 
thereby improving aggregate stability in a field experiment at Central China (Dai et al., 2015). 
Meanwhile, sustainable farming management including no till age, continuous crop coverage and 
light pruning increased the content of water stable aggregates by enhancing the accumulation of 
glomalin in macro-aggregates in an olive grove at southern Italy (Lombardo et al., 2019). Binding 
agents in macro-aggregates contributed more to soil structure than those in other sized aggregates, 
due to its higher content and sensitivity to environmental changes (Singh et al., 2016). 

The content of GRSP and SOC in aggregates is directly influenced by soil properties and plant 
characteristics (Clair and Lynch, 2010; Wang et al., 2016; Wang et al., 2018; Yue et al., 2018). For 
example, soils with neutral or low pH improve root and fungal growth, and produce abundant 
glomalin and SOC (Wang et al., 2015; Zhong et al., 2017). Low content of soil nutrients, such as 
carbon (C), nitrogen (N), and phosphorus (P) increases the lifespan of ectomycorrhizal roots and 
reduces the rate of GRSP decomposition (Lovelock et al., 2003, 2004). Increased plant biomass is 
beneficial to arbuscular mycorrhizal fungi (AMF) colonization and GRSP accumulation (Wang et 
al., 2018). Climatic conditions, such as precipitation and temperature, have profound effects on 
soil properties and plant characteristics, which may lead to changes in GRSP and SOC (Treseder 
and Allen, 2000). Therefore, climatic conditions change soil structure by affecting these factors, 
but the mechanisms remain unclear. On one hand, the increasing precipitation and temperature 
across semi-arid and semi-humid forest ecosystems improves soil quality, which, in turn 
accelerates the decomposition of organic matter (Clair and Lynch, 2010). On the other hand, plant 
biomass and primary productivity are higher in semi-humid forests than in semi-arid forests, 
which are conducive to AMF infection and glomalin production (Wang et al., 2016; Yue et al., 
2018). The decomposition and production of binding agents are enhanced at the same time, and 
soil aggregate stability depends on the combined effects of these biochemical processes along a 
climatic gradient. Based on previous studies, this work will explore the variations of binding 
agents and soil aggregate stability caused by climate change. The effects of plant characteristics 
and soil properties during this process were evaluated. 

In this study, Robinia pseudoacacia L. forests along a climatic gradient (from semi-arid to 
semi-humid) on the Loess Plateau, China, were selected as the research objects. We focused on 
the biomass of R. pseudoacacia, soil physical-chemical properties and distribution of EE-GRSP, 
total GRSP (T-GRSP) and SOC in soil aggregates. Using correlation analysis and combining with 
redundancy analysis (RDA), we analyzed the contribution of binding agents to soil aggregate 
stability and the effects of plant biomass and soil physical-chemical properties on the stability. 


2 Materials and methods 


2.1 Study area 


R. pseudoacacia forests at Shenmu City (SM), Suide County (SD), Ansai District (AS) and 
Chunhua County (CH) from the north to the south of the Loess Plateau, China were selected 
(34°46'-38°44'N, 108°31'-110°29’E; 1100-1200 m a.s.l.). Sampling locations are presented in 
Figure 1. Annual average temperature of the four study sites was 9.2°C—10.6°C, annual minimum 
temperature was —24.0°C— -—29.0°C in January, and annual maximum temperature was 
36.0°C-42.0°C in July. Average annual precipitation was 390.3-585.7 mm. The study sites 
belonged to the transitional zone from semi-arid to semi-humid forest ecosystems and was 
dominated by continental monsoon climate, with distinct seasonal and micro-climates. The soil of 
R. pseudoacacia forests was mainly loessal soil with a slope of 25°—30°. The canopy density of 
the forests was 0.5—0.7 and the diameter at breast height (DBH) was 10-20 cm. The shrubs and 
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herbs in the forests included Ziziphus jujuba, Rosa xanthina, Stipa bungeana, Carex lanceolata, 
Artemisia scoparia, Saussurea romuleifolia. Information of the four study sites is shown in Table 1. 
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Fig. 1 Location of the study sites. SM, Shenmu City; SD, Suide County; AS, Ansai District; CH, Chunhua 
County. 


Table 1 Information of the study sites 


Tree 


Site Latitude and = Altitude Precipitation Temperature Evaporation Main height DBH Soil 
longitude (m) (mm) (°C) (mm) plants ts k (cm) type 
Locust 

38°44'46"N j Sand 

SM 110°29'48"E 1204.4 390.3 9.2 1336.6 s 6.91 18.18 loess 
37°25'19"N Locust, ma 

SD 110°09'55"E 1178.3 410.5 10.1 2061.9 Miscanthus 7.44 16.33 Loess 
36°44'25"N Locust, ; 

AS 109°15’11"E 1110.9 506.5 9.2 1000.0 Miscanthus 12.35 17.35 Loess 
34°46'48"N Locust, : 

CH 108°31'44"E 1276.5 585.7 10.6 520.0 Jujube 9.28 15.93 Loess 


Note: SM, Shenmu City; SD, Suide County; AS, Ansai District; CH, Chunhua County. DBH, diameter at breast height. 
2.2 Experimental design 


In September 2017, we selected plots (20 mx10 m in each plot) with the same direction and slope 
on the basis of site conditions. Each site of R. pseudoacacia forests had five plots that were 
separated at least 20 m apart from each other. Species composition, height and DBH of R. 
pseudoacacia in each plot were recorded. Three undisturbed soil samples were randomly 
collected from 0-20 cm soil layer after litter was removed. The soil samples were frozen and 
carefully transported back to the laboratory. The three soil samples in each plot were mixed and 
gently broken apart. The roots and stones were removed simultaneously. One part of the soil 
sample was frozen for determining water content and AMF infection rate. Another part was used 
for determining soil physical-chemical properties, aggregate composition, SOC and GRSP 
contents after air drying. 


2.3 Soil aggregate distribution 
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Soil aggregate distribution was determined using the wet sieving method (Six et al., 1998). Filter 
paper was spread on a glass dish and soaked with deionized water. Approximately 100 g soil 
sample was evenly spread on the filter paper for 30 min. The soil sample was submerged on the 
water surface, and the impurities on the water surface were removed. Then, the subsample was 
shaken for 1 min on a wet sieve machine with a frequency of 50 times/min and amplitude of 3 cm. 
The subsample was separated into three classes for the different sized aggregates: >0.250 mm 
(macro-aggregate), 0.250—0.053 mm (micro-aggregate) and <0.053 mm (fine material). All 
separated soil and water were oven-dried at 45°C, weighed and saved for laboratory analysis. 


2.4 GRSP extraction and determination 


GRSP content was determined in accordance with the method of Wrightand Upadhyaya (1998). 
Soil sample (1 g) and 8 mL of 50 mmol/L citrate solution (pH=7.0) were added to a 10-mL 
centrifuge tube and autoclaved at 121°C for 60 min. The sample was centrifuged at 10,000 r/min 
for 5 min, and the supernatant was poured into a 50-mL centrifuge tube. This procedure was 
repeated five times on the same soil sample to obtain the T-GRSP extract. The extract was mixed 
and stored for 3 d at 4°C. EE-GRSP extract was obtained in accordance with the above mentioned 
method, but the procedure was repeated once. T-GRSP and EE-GRSP contents of the extract were 
tested spectrophotometrically at 590 nm light via the Bradford dye-binding assay with bovine 
serum albumin as the standard. 


2.5 Soil property determination 


Dichromate oxidation (Mebius, 1960) and Kjeldahl methods (Page et al., 1982) were used to 
measure SOC and N contents, respectively. Soil P content was measured by persulfate oxidation 
followed by colorimetric analysis (Lu, 2000). Soil pH was determined by potentiometry (Lu, 
2000). The infection rate of AMF was observed via the trypan blue coloration method (Phillips 
and Hayman, 1970). Soil water content was determined by the drying method (Lu, 2000). 


2.6 DBH and biomass determination of R. pseudoacacia 


Above-ground biomass, root biomass and total biomass were determined using the allometric 
growth model of R. pseudoacacia forests on the Loess Plateau, China (Li, 2015). DBH was 
measured using a growth cone to extract the pith, and dichotomy was used to estimate the height 
of trees. 
The allometric growth model was as follows: 

In(M)=a+bIn(DBH), (1) 
where M is the biomass (kg); a and b are the estimates of model parameters that vary according to 
organ and diameter of trees, respectively; and DBH is the diameter at breast height (cm). 


2.7 Statistical analyses 


Soil aggregate stability was evaluated by mean weight diameter (MWD) and geometric mean 
diameter (GMD). MWD and GMD were calculated using the following equations: 


MWD=)""_ XW; > (2) 


GMD=exp ad , (3) 
a 

where x; is the mean diameter of each sized aggregate (mm); w; is the percentage of each sized 

aggregate (%); and n is the number of different sized aggregates. 

One-way variance analysis was used to analyze the differences in soil physical-chemical 
property, aggregate distribution, stability, and SOC and GRSP contents among four sites (P<0.05). 
Correlation analysis was used to compare the relationship between aggregate stability and binding 
agents. The above mentioned analysis was conducted using SPSS 20.0 statistical software 
package (SPSS Inc., Chicago, IL, USA). RDA was used to evaluate the effects of soil 


physical-chemical property and plant biomass on binding agents in soil aggregates. Analysis was 
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conducted by CANOCO (v4.5, Microcomputer Power, Ithaca, NY, USA). 


3 Results 


3.1 Soil property and R. pseudoacacia biomass 


N, SOC and water contents were significantly increased, and pH was significantly decreased 
along the climatic gradient (P<0.05; Table 2). pH at the CH site was almost neutral. AMF 
infection rate initially decreased and then increased along the climatic gradient. SM and CH sites 
had the highest value of infection. Above-ground biomass, root biomass and total biomass of R. 
pseudoacacia significantly increased along the climatic gradient and had the highest value at the 
CH site. 


Table 2 Soil property and R. pseudoacacia biomass of the study sites 


idek Study site 
SM SD AS CH 

N (mg/g) 0.310.059” 0.320.057” 0.330.033” 0.6940.325° 

P (mg/g) 0.510.017” 0.550.006% 0.57+0.030* 0.56+0.050° 

SOC (mg/g) 4.73+0.759° 3.01+0.349° 4.62+0.868° 9.63+4.791" 

pH 8.40+0.034° 8.48+0.039° 8.41+0.048° 8.21+0.052° 

Water content (%) 6.36+0.736° 10.20+1.319° 10.7541.114° 16.53+1.949* 

AMF infection rate (%) 83.94+4.698" 35.6144.525° 42.63+12.401° 77 A0#4.565° 
Above-ground biomass (kg) 76.8445.396° 60.11+20.442° 102.73+421.453* 111.63423.854* 

Root biomass (kg) 39.03+42.488° 31.114+9.675° 63.76+11.670° 52.95+9.772* 
Total biomass (kg) 115.8747.884° 91.22430.116° 166.49433.122* 164.58+33.626° 


Note: SM, Shenmu City; SD, Suide County; AS, Ansai District; CH, Chunhua County; N, nitrogen; P, phosphorus; SOC, soil organic 
carbon; AMF, arbuscular mycorrhizal fungi. Different lowercase letters within the same row indicate significant differences among the 
four study sites at P<0.05 level. Meant+SD. 


3.2 Distribution and stability of soil aggregates 


The proportions of macro-aggregates, micro-aggregates and fine materials were in the range of 
19.43%—-710.44%, 16.12%-47.61% and 4.34%—30.44%, respectively. Soil at the CH site was 
dominated by macro-aggregates, whereas soils at the SM, SD and AS sites were primarily 
micro-aggregates (Fig. 2). From semi-arid to semi-humid forests, the proportion of 
macro-aggregates significantly increased, whereas those of micro-aggregates and fine materials 
significantly decreased (P<0.05). MWD and GMD were significantly enhanced along the climatic 
gradient (P<0.05). 


3.3 Distribution of SOC and GRSP in soil aggregates 


From semi-arid to semi-humid forests, the contents of T-GRSP and SOC in macro- and 
micro-aggregates significantly increased and had the highest values at the CH site (P<0.05; Fig. 
3); the contents of T-GRSP and SOC in fine materials significantly decreased and had the highest 
value at the SM site; the content of EE-GRSP in macro- and micro-aggregates initially decreased 
and then increased with the highest values at the SM and CH sites, whereas the content of 
EE-GRSP in fine materials significantly decreased. Across all sized aggregates, the content of 
GRSP in macro-aggregates had the highest value, followed by that in micro-aggregates and fine 
materials. The distribution of SOC in aggregates was consistent with that in GRSP, except for that 
at the SM site. 
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Fig. 2 Distribution (a) and stability (b) of soil aggregates at the study sites. Different lowercase letters indicate 
the significant differences among the study sites for each parameter at P<0.05 level. Different uppercase letters 
indicate significant difference among different sized aggregates at P<0.05 level. SM, Shenmu City; SD, Suide 
County; AS, Ansai District; CH, Chunhua County; MWD, mean weight diameter; GMD, geometric mean 
diameter; macro-aggregates, >0.250 mm; micro-aggregates, 0.250—0.053 mm; fine materials, <0.053 mm. 
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Fig. 3 Distribution of easily extractable glomalin-related soil protein (EE-GRSP, a), total extractable 
glomalin-related soil protein (T-GRSP, b) and soil organic carbon (SOC, c) in aggregates. Different lowercase 
letters indicate significant differences among the study sites for each parameter at P<0.05 level. Different 
uppercase letters indicate significant differences among different sized aggregates at P<0.05 level. SM, Shenmu 
City; SD, Suide County; AS, Ansai District; CH, Chunhua County; macro-aggregates, >0.250 mm; 
micro-aggregates, 0.250—0.053 mm; fine materials, <0.053 mm. 


3.4 Correlation between soil aggregate stability and binding agents 
GRSP and SOC in macro-aggregates and micro-aggregates were significantly related with MWD 
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(P<0.05, Table 3), whereas no significant relationship was observed between MWD and binding 
agents in fine materials (P>0.05). The correlation coefficient between SOC and MWD was 
greater than that between T-GRSP and MWD and that between EE-GRSP and MWD. The 
correlation coefficient between MWD and binding agents in macro-aggregates was greater than 
that in micro-aggregates. 


Table 3 Bivariate correlation between mean weight diameter (MWD) and binding agents in soil aggregates 


Index Bulk soil Macro-aggregates Micro-aggregates Fine materials 
EE-GRSP ND 0.574" 0.339 -0.001 
T-GRSP ND 0.711” 0.587" 0.110 
soc 0.750" 0.795” 0.758” 0.115 


Note: EE-GRSP and T-GRSP are easily extractable and total glomalin-related soil proteins, respectively; SOC, soil organic carbon; 
macro-aggregates, >0.250 mm; micro-aggregates, 0.250-0.053 mm; fine materials, <0.053 mm; ND, not determined. “and “ indicate 
significances at P<0.05 and P<0.01 levels, respectively. 


3.5 RDA of ecological factors and binding agents 


The first and second axes accounted for 96.65% of the total variance in the data set (Fig. 4). The 
effect of soil property on binding agents was greater than that of R. pseudoacacia biomass. AMF 
infection rate had a strong effect on T-GRSP in macro-aggregates. SOC in macro-aggregates was 
affected by N, water content and total biomass. pH had negative effects on T-GRSP and SOC in 
macro-aggregates. 
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Fig. 4 Redundancy analysis (RDA) of binding agents, soil property and R. pseudoacacia growth index. N, 

nitrogen; P, phosphorus; SOC, soil organic carbon; AMF, arbuscular mycorrhizal fungi; EE-GRSP and T-GRSP 

are easily extractable and total glomalin-related soil proteins, respectively; macro-aggregates, >0.250 mm; 

micro-aggregates, 0.250—0.053 mm; fine materials, <0.053 mm. 


4 Discussion 


4.1 Composition and stability of soil aggregates along the climatic gradient 


Soil macro-aggregate proportion increased from semi-arid to semi-humid R. pseudoacacia forests 
on the Loess Plateau, China. This result was due to the increased content of binding agents and 
suitable precipitation and temperature conditions. Previous studies have confirmed that climatic 
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condition promoted the large aggregate formation and binding agent production while decreasing 
the micro-aggregate proportion (Zhang et al., 2014; Xie et al., 2015; Luna et al., 2016). Different 
sized aggregates had varied changes because the increased soil quality caused by climate 
improvement converted small sized aggregates into large sized aggregates, thereby decreasing the 
proportion of micro-aggregates. Dai et al. (2015) also found that the proportion of 2-50 um 
aggregates and binding agents had different responses to experimental treatments. In this study, 
soil aggregate stability was enhanced from semi-arid to semi-humid forests, which was related to 
the various effects of binding agents on soil aggregates. 


4.2 Distribution of binding agents in soil aggregates 


Improvement of climatic conditions (such as increasing precipitation and temperature) enhanced 
GRSP production and promoted the conversion of GRSP from fine materials to micro-aggregates 
and macro-aggregates. This phenomenon was consistent with the aggregate formation model, 
which was proposed by Six et al. (2004). That is, large sized aggregates were formed by the 
combination of small sized aggregates and binding agents. These results indicated that the size of 
aggregates should be considered when climatic conditions affected GRSP in forest ecosystems on 
the Loess Plateau, China. Different from farmland, GRSP in all sized aggregates of R. 
pseudoacacia forests had a similar response along the climatic gradient, which might be related to 
the infertile soil, limited precipitation and low temperature (Zhang et al., 2014). In addition, our 
results showed that the content of GRSP at the SM site (semi-arid) was significantly higher than 
those at the SD and AS sites. First, this phenomenon was attributed to the relatively poor soil at 
the SM site, which was conducive to the AMF growth, prolonging the lifespan and decreasing the 
decomposition rate of organic matter (Lovelock et al., 2003; Lovelock et al., 2004). Second, 
precipitation and soil water content are low at the SM site, therefore binding agents, such as 
GRSP, could not further combine with soil particles to form aggregates, resulting in no 
improvement in aggregate stability. As shown in Figure 2, MWD and macro-aggregate 
proportions at the SM site were significantly lower than those at other sites. These results indicate 
that aggregate formation by glomalin requires suitable hydrothermal conditions on the 
semi-arid/semi-humid forests. 

SOC and GRSP in aggregates had different responses along the climatic gradient. Low annual 
precipitation and temperature resulted in poor soil quality and less production of SOC, whereas 
climatic condition at the SM site was favored for GRSP accumulation. Meanwhile, GRSP had 
lesser contribution to SOC compared with other components (Zhu et al., 2019). A study in 
Northeast China had found that SOC was primarily composed of inert carbon during urbanization, 
and glomalin did not contribute to the accumulation of inert carbon (Wang et al., 2018). In this 
study, SOC in aggregates was primarily affected by climatic conditions, whereas GRSP was 
affected by climate conditions and AMF characteristics. 


4.3 Relationship between binding agents and aggregate stability 


Binding agents in macro-aggregate had more contribution to aggregate stability than those in 
micro-aggregate and fine materials. From semi-arid to semi-humid forests, climatic conditions 
were favored for combination of macro-aggregate and improving soil structure (Zhang et al., 
2014). 

The contribution of SOC to aggregate stability was greater than those of T-GRSP and EE-GRSP, 
which was not in line with our assumption. The effect of binding agents on aggregate stability 
depended on two aspects: aggregation efficiency and content of binding agents. These two aspects 
determined the relationship between binding agents and soil stability (Rillig et al., 2001; Six et al., 
2002; Rillig, 2004). EE-GRSP had a high binding efficiency, but considerably low content, which 
resulted in negligible influence on aggregate stability. Zhu et al. (2019) did not find a significant 
relationship between EE-GRSP and aggregate stability in southwestern China. SOC contained 
various inert components with a low binding ability, but its content was higher than that of 
EE-GRSP (SOC content in this study was 33.8 times higher than that of EE-GRSP). Therefore, 
SOC had a greater influence on aggregate stability than EE-GRSP. 
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4.4 RDA of ecological factors and binding agents 


From semi-arid to semi-humid forests on the Loess Plateau, climatic conditions promoted the 
growth of R. pseudoacacia and soil quality (Clair and Lynch, 2010; Wang et al., 2016; Yue et al., 
2018). AMF infection rate was high at the SM site (semi-arid) because soil with low nutrients was 
conducive to fungal growth and prolonged mycorrhizal lifespan (Treseder and Allen, 2000). 

RDA showed that the effect of soil property on binding agent was greater than that of plant 
biomass. This result was consistent with the research in a farmland (Sarker et al., 2018). Although 
plant biomass of R. pseudoacacia forests was evidently superior to that of farmland, soil property 
was still the most important factor affecting binding agent in soil aggregates. Soil N was the 
important index that increased plant growth and microbial activity, leading to a strong effect on 
the distribution of SOC in aggregates (Aber et al., 1998; Hong et al., 2017). Value of pH 
determined the rate of decomposition of litter and activity of AMF. Neutral soil was beneficial to 
the accumulation of different binding agents (Wang et al., 2018). Infection rate directly controlled 
AMF colonization and GRSP production. Soil water content was a limiting factor for microbial 
community in the sites of semi-humid and semi-arid forests, resulting in a remarkable effect on 
the accumulation of OC in aggregates. Soil N, pH, AMF infection rate and water content had 
positive effects on accumulation of binding agents. The abovementioned results were supported 
by previous studies. For example, Wang et al. (2014) showed that soil pH was the primary factor 
that changed glomalin between soil layers. C, N and water content were important for GRSP 
accumulation in different soil types of farmlands (Zhong et al., 2017). From semi-arid to 
semi-humid forests, our study emphasized the role of soil property in improving stability of soil 
structure. 


5 Conclusions 


Plant primary production, soil chemistry and structure in R. pseudoacacia forests gradually 
changed along the climatic gradient on the Loess Plateau, China. This study tested the mechanism 
of soil structure changing from semi-arid to semi-humid. First, SOC was a more efficient binding 
agent than that previously suggested. Second, soil chemistry had greater contribution to soil 
structure than vegetation factor. Third, some indices, such as N and AMF infection rate were 
emphasized during aggregate formation, providing reference for future research. These results 
revealed a potential way through which environmental conditions affected aggregate stability and 
multiple parts therein were elucidated. 
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